High efficiency transfer of quantum information and multi-particle entanglement generation in 

translation invariant quantum chains. 
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We demonstrate that a translation invariant chain of interacting quantum systems can be used for high effi- 
ciency transfer of quantum entanglement and the generation of multi-particle entanglement over large distances 
and between arbitrary sites without the requirement of precise spatial or temporal control. The scheme is largely 
insensitive to disorder and random coupling strengths in the chain. We discuss harmonic oscillator systems both 
in the case of arbitrary Gaussian states and in situations when at most one excitation is in the system. The latter 
case which we prove to be equivalent to an xy-spin chain may be used to generate genuine multi particle entan- 
glement. Such a 'quantum data bus' may prove useful in future solid state architectures for quantum information 
processing. 

PACS numbers: 03.67.-a,03.67.Hk 



The realization of quantum communication and computa- 
tion requires at various stages the mapping between station- 
ary and flying qubits and the subsequent transfer of quantum 
information between different units of our quantum informa- 
tion processing devices. Traditionally the stationary forms of 
qubits are massive systems such as atoms, ions, quantum dots 
or Josephson junctions while the flying qubit is a photon, ie 
radiation. Photons might be optimal when considering long 
distance communication where they may travel through free 
space or optical fibres. In very small quantum information 
processing devices such as condensed matter systems how- 
ever, this is difficult as the length scale both of the compo- 
nent parts and their separation will generally be below op- 
tical wavelengths. In this situation, it is worth considering 
novel approaches for the communication of quantum infor- 
mation and the generation of entanglement. To this end it 
is of interest to consider the properties of interacting quan- 
tum systems and here in particular those of harmonic sys- 
tems that are realized in various condensed matter physics 
settings such as nano-mechanical oscillators. While static har- 
monic (or spin) systems near their ground state do not exhibit 
long distance entanglement the situation changes drasti- 
cally when consideringtime-dependent properties of interact- 
ing quantum systems Indeed, solid state devices such as 
arrays of nano-mechanical oscillators, described as interacting 
harmonic oscillators, allow for the generation transfer and 
manipulation of entanglement 01 with a minimum of spatial 
and temporal control. However, in translation invariant sys- 
tems the efficiency for this transfer decreased with distance. 
This can be overcome either by making the coupling strengths 
between neighboring systems position dependent 00] or by 
active steps such as quantum repeater stages |0] or conclu- 
sive transfer |j7|. Nevertheless, active steps or the fabrication 
of precisely manufactured spatially dependent couplings are 
difficult in practice and will require a significant degree of 
control. Furthermore, the precise value of the coupling pa- 
rameters and the timing of the operations will depend on the 
distance across which one aims to transfer quantum informa- 



tion. Consequently, it would be desirable to achieve high ef- 
ficiency transmission of quantum information between arbi- 
trary places and distances with minimal spatial and temporal 
control. In the following we show that this is indeed possible 
employing translation invariant chains of interacting quantum 
systems with stationary couplings. 

We first describe the system, termed quantum data bus, and 
demonstrate its functionality by numerical examples. Then 
we present an approximate analytical model that reveals the 
basic physical mechanism that is utilized in the operation of 
the quantum data bus. This model then allows us maximize 
entanglement transfer efficiency and transmission speed of the 
quantum data bus by adjusting the eigenfrequencies of the 
sender and receiver system. It also explains why the trans- 
mission is largely insensitive to disorder and random coupling 
strengths in the ring. We discuss the scaling behavior of the 
time that is required for the transfer between distant sites at a 
given efficiency. Finally, we show that in the regime where at 
most one excitation is in the system, a quantum data bus made 
of interacting harmonic oscillators becomes equivalent to an 
interacting spin chain. A further application for such a chain 
is the generation of three or multi-particle entangled states as 
we also show in this paper. The model we have in mind is de- 
picted in Fig. 1. A ring of interacting quantum systems (blue 
circles) forms the quantum data bus. At arbitrary positions 
on the ring two further quantum systems (red) may be cou- 
pled weakly to the ring. The subsequent time evolution will 
allow for high efficiency transfer of entanglement between the 
two distinguished quantum systems. The following ideas are 
not restricted to the specific ring like geometry presented here. 
The depicted geometry has been chosen because it simplifies 
the analytical treatment of the problem. 

So far we have not specified any particular type of quantum 
systems nor their mutual interaction. In the following we will 
consider as an example the case of coupled harmonic oscilla- 
tors which may be realized in various solid state settings 1 3] . 
Towards the end of this article we will also consider a setting 
which is equivalent to a spin chain with an xy-interaction. Let 
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FIG. 1 : A ring of interacting quantum systems (blue circles) forms 
the quantum data bus. Two further quantum systems a and b (red) 
may couple weakly at arbitrary positions to the ring. The subsequent 
time evolution will allow for the transfer of quantum information 
or entanglement that exists between systems a and the decoupled 
system c from system a to system b. System c is decoupled but might 
be entangled with a to study transfer of entanglement. 



us assume that the ring consists of M harmonically coupled 
identical oscillators. We set h = oj ~ m = 1 and denote the 
coupling strength in the ring by c. We assume that oscilla- 
tor a (b) couples to oscillator 1 (k) with strength e so that the 
Hamilton operator of this system including the oscillators a, b 
and c is given by 
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with the potential matrix V given by Vkk ~ 1 + 2c and 
Vk^k+i = Vk+i,k = and Vim = Vmi = -c for all k 
and zero otherwise. 

In the following we will demonstrate that it is indeed pos- 
sible to transmit quantum information through this system 
with high efficiency but minimal spatial and temporal con- 
trol. For the moment we are focussing on Gaussian states, 
i.e., states whose characteristic function or Wigner function 
is Gaussian. The characteristic function determines the quan- 
tum state and as any Gaussian is determined by its first an 
second moments the same applies to the corresponding quan- 
tum state (si]. In the present setting the first moments will 
not be directly relevant (they correspond to biases that can 
be removed by redefining the coordinate origin). Therefore 
the state of the system is determined by second moments that 
can be arranged in the symmetric 2N x 2Af-covariance ma- 
trix Tr^s = 2Re((i? - {R)){S - {S))), where R and S 
stand for the canonical operators xi, . . . ,Xn and pi, . . . ,pn- 
Employing the Hamilton operator eq. Q we can now study 
numerically the quality of the entanglement transfer between 
oscillators a and b. Let us consider the situation where the 
harmonic oscillator a and c are initially in a pure entangled 
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FIG. 2: The efticiency of entanglement transfer, defined as ratio of 
transmitted entanglement to initial entanglement, between the oscil- 
lator a and b at opposite ends of a quantum data bus consisting of 20 
oscillators and a nearest neighbor coupling strength of c = 1. The 
coupling strength of the oscillators a and b to the quantum data bus 
is e = 0.015 (dashed line) and e — 0.021 (solid line). The speed 
of propagation is proportional to the coupling strength e and in both 
cases we observe maximal efficiency exceeding 0.99. 



two-mode squeezed state 



(2) 



with q = tanh | which possesses the covariance matrix 
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The entanglement as quantified by the logarithmic negativity 
of this state is then given by £'Ar(|-0)) = 2r. The time 
evolution of the entanglement between the oscillator a and b 
at opposite ends of a quantum data bus consisting of 20 os- 
cillators and a nearest neighbor coupling strength of c = 1 is 
given in Fig. |2] The propagation speed is proportional to e and 
the efficiency of entanglement transfer decreases weakly with 
increasing speed. In both cases we find a maximal efficiency, 
defined as the ratio of transmitted entanglement to initial en- 
tanglement, exceeding 0.99. It should be noted that the time 
required for the transfer of entanglement is independent of the 
distance of places where the oscillators a and b couple to the 
ring. The high quality of the entanglement transfer and its 
independence on the position of sender and receiver will be 
successfully explained by the following model that encapsu- 
lates the essential physics in the system. 

We first observe that the unitary matrix fl with matrix ele- 
ments 
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achieves V = il^A^fl with a diagonal matrix A such that 
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1 + 2c — 2c cos- 
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Then we can define the normal mode variables 
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which ensure that [X^, Pi] = [xk,pi] = iSki, ie the canonical 
commutation relations are preserved. Note that we will use 
the convention Xq = Xm and Pq = Pm which reflects the 
periodic boundary conditions of the quantum data bus. Fur- 
ther denote Xa^b,c = Xa,b,c and Pa,b,c = Pa,b,c to make the 
notation more uniform. In these normal modes we can write 
the Hamiltonian eq. Q as 
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Defining the annihilation operators 



Ak 



KkkXk + iPl 



V2A 
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kk 
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we can rewrite the Hamilton operator in terms of the Ak ■ In- 
deed, shifting the zero of energy and moving to an interaction 
picture with respect to 



1 1 
H^:^J2PlPk + AlkXlXk + - P^ + ^f 
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we find 
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Hj = AiAa + -AlAb~eJ2 
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In this interaction picture we find for e = that Al{t) = {X^, Xa, Xb, Xm)'^ we find 
Alio) while for finite e we have 4^i(t) = i[Hi, Al{t)]. In 



time dependent perturbation theory using 



Ai{t + At)-Ai{t) 
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ds,[Hj{s,),Ai{t)] 



and setting i ^ At <C ^ we find that to first order in e 
the modes described by Aa and Ab are only coupled to one 
collective mode, namely the center of mass mode, described 
by Am- Therefore we can ignore the contributions from all 
other eigenmodes. Shifting the zero of energy again we finally 
obtain the simplified set of equations 



H. 



approx 



\AiAa + ^AlAb + ^AjAm 
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[iAa + Ab)Al + {Ai + Al)AM].{8) 



Now we write this Hamiltonian again in the quadrature com- 
ponents. Defining P = {Pc, Pa, Pb, Pm)'^ and X = 



where 
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The above set of approximate equations of motion can be un- 
derstood intuitively by a simple mechanical model. Indeed 
they describe the motion of a very heavy central pendulum 
(corresponding to the oscillators in the quantum data bus) 
that is coupled weakly to two comparatively light oscillators. 
From undergraduate mechanics we know that if one of the 
light pendulum is originally oscillating then, after some time, 
it will have stopped oscillating, while the other light pendu- 
lum is now oscillating with almost the same amplitude while 
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the heavy central pendulum remains essentially at rest. From 
this simple mechanical picture the dynamics in the quantum 
setting that is described below can be understood quite intu- 
itively. 

That the above approximate Hamiltonian represents a good 
approximation to the true dynamics can be seen from a com- 
parison of the exact time evolution with that generated by the 
approximate Hamiltonian. In figure |3] we chose a ring with 
20 oscillators, c = 10 and e = 0.021. The observable mis- 
match of between the frequencies is due to second order cor- 
rections to the approximate model. Furthermore, as the exci- 
tation of a quantum in the center of mass mode corresponds 
to the simultaneous in-phase motion of all the oscillators in 
the quantum data bus and the fact that the oscillators a and 
b couple predominantly to the center of mass mode M also 
explains why the entanglement transfer is distance indepen- 
dent. The frequency of the center-of-mass mode is indepen- 
dent of both the coupling strengths between oscillators in the 
chain and possible disorder in which oscillators are coupled 
as long as there is a connected path of coupled oscillators be- 
tween a and b. Therefore the transmission of quantum infor- 
mation between a and h will depend only weakly on disorder 
and randomness. The main correction arises when the fre- 
quency separation between the lowest two eigen-modes be- 
comes small so that off-resonant couplings of the oscillators 
a and b to modes other than the center-of-mass mode become 
non-negligible. In the derivation of eq. (|8j we have neglected 
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FIG. 3: Comparison of the approximation and the exact time evo- 
lution for a ring of 20 oscillators with c = 10 and e = 0.021. The 
dotted line is the approximation and the solid line is exact. 

many terms that led to oscillating contributions in the Hamil- 
tonian. These neglected terms will lead to corrections whose 
size will depend on the length M of the quantum data bus 
and will affect both, the speed of propagation but also its ef- 
ficiency. As we can always adjust the waiting time, correc- 
tions to the propagation speed are less relevant. An efficiency 
reduction due to population losses is more serious as it will 
require the application of error correction methods. For the 
setting described by the Hamiltonian eq. Q we will now 



present an estimate the size of these corrections. Indeed, fol- 
lowing eq. we neglect all terms that couple the modes a 
and 6 non-resonantly to eigenmodes different than the center 
of mass mode. For small coupling strength e the rapid oscilla- 
tions will reduce the population in these modes significantly. 
The mode, other than the center of mass mode, with the small- 
est oscillation frequency will be the mode M—1. For large M 
the frequency difference to the center of mass mode is given 
by A = 27r^c/M^ and the coupling strength to this mode is of 
the order of e/2\/M . The loss of population into this mode, 

which is of order ( ) | ] , should be small, ie 



e 47r-^ 
c A/3/2 



(9) 



From eq. (|8|i we determine that the speed of propagation is 
proportional to e/ (2Af ). For a prescribed transfer efficiency 
the transmission time therefore scales as 



T . 



2M 
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(10) 



The error source that enters this scaling is the population loss 
to off-resonant modes. This suggest that this scaling can be 
improved considerably when one allows for a fixed frequency 
difference of the oscillators a and b compared to the oscilla- 
tors in the quantum data bus such that the oscillators a and b 
become resonant with a different collective mode whose fre- 
quency difference to the neighboring modes is as large as pos- 
sible. Indeed, if we couple a and 6 to the mode described by 
Am/^, ie if we shift the eigenfrequencies of a and b by 2c, 
then we find that the next mode is separated by a frequency 
difference A = 2c-k/M and the couplings strength to this 
neighbouring mode is again of the order of e/2\/M. As a 
consequence the population loss to other modes is of the or- 
der of (e-y/M / (47rc))^ and we only need to ensure that 



e 47r 
- < ^= 

c 



(11) 



so that the propagation time for fixed transfer efficiency scales 
as 



2M 
e 



~2ttc 



(12) 



This improved scaling has been achieved by coupling the os- 
cillators a and b, but it should be noted that unlike the center of 
mass mode M this mode will have nodes. As a consequence 
there will be relative positions of oscillator a and b such that 
they will not couple, namely, when one is sitting at a node 
while the other is at an anti-node. Indeed, the A// 4 mode has 
a node at every second oscillator of the ring, so that the os- 
cillators a and b couple only when their distance is an even 
number of oscillators. 

Finally we demonstrate that the above considerations are 
not restricted to the continuous variable regime and the prop- 
erties of Gaussian continuous variable states. Indeed, we can 
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equally well consider a situation in which we restrict our dy- 
namics to the subspace spanned by the the vacuum and those 
states that corresponds to a single excitation. The deriva- 
tion of the approximate model starting from Hamiltonian eq. 
([0 remains of course valid. In the basis represented by 

the states {|00)a6|0)A/, |10)a6|0)M, |01)a6|0)M, |00)a6|l)M} 

where \l)]\f = X^rli '^J-l'^'^'^""™) rewrite the 

Hamiltonian eq. (|8jl as 



Ha 



■pprox 



2\M 2 ' ^ 



(13) 



Ml 



which corresponds to a spin chain with an xy-interaction fioll 
under the same constraint of considering at most a single ex- 
citation. This similarity between a harmonic oscillator sys- 
tems and a spin chain is not due to the approximations in the 
derivation of eq. (|8j but a generic feature when one limits the 
number of excitations to at most one. A simple computation 
shows that generally a harmonic chain in the rotating wave- 
approximation and the single excitation regime will be equiv- 
alent to a spin chain with xy-interaction in the same regime. 

Continuing in this setting of a spin chain we now demon- 
strate that multi-particle entanglement 1 11] can be generated 
with a quantum data bus extending the ideas employed in the 
paper so far While we focus on the discrete case, i.e. the spin 
Hamiltonian, one may carry out a similar investigation for the 
harmonic oscillator case. For the purpose of the generation of 
entangled states there is no need of including the decoupled 
oscillator in the discussion. While the ideas presented below 
may easily be generalized to many oscillators let us, for sim- 
plicity, consider three oscillators coupled to the chain. We 
assume that the oscillator a {b, c) couples to oscillator 1 (fci, 
of the quantum data bus with strength e so that the Hamil- 
ton operator of this system is a generalization of eq. (I13> and 
is given by 



3-M 
M 



, as well as 



'12 + ^2 
C{M,t) = cos(\/^^r/2) 



/l2 + (1-2 

S{M,r) = sin(^^^r/2) 
the time dependent coefficients are given by 

-iar/2 



(17) 



a(r) = i 
6(r) 



Vl2 + Ma2 

— iOLT jl 
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Vl2 + Ma^ 
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c(t) - d(r) = -- + — 



C{M, t) + ^ 



S{M, 



Vl2 + Ma^ 



From eq. (I18> we can see that in the limit of large M, the 
coefficient a(T) tends to zero indicating that the quantum data 
bus disentangles from the three oscillators and the result may 
be a W state of the form 

\W,,y,,)®\0)M - (x|100)+y|010) + z|001))®|0)M (19) 

In order to see that this is the case we plot the overlap between 
the W state eq. il9i and the evolved state obtained with the 
use of eq. (I16> and eq. ( I18> . For the case of 70 oscillators in 
the chain, and the W state with x = —y = —z = 1/ \/3, it is 
shown in figure|4]that the overlap is about 96%. The overlap is 
not complete due to fact that the population of the state |0001) 
does not vanish and that the coefficients eq. ( I18> have a small 
imaginary part. It should be noted that the second of these 
effects can be reduced significantly by employing subsequent 
local unitary rotations which are irrelevant when considering 
entanglement properties of a state. In figure |5l we show the 
populations in function of the scaled time r. For generating 



(14) 



Let us now suppose that the initial state of the system is 



Ml 



0.9- 



0.8- 



> 

O 0.7- 



|^A(0)) = |l)a|0)b|0)c|0)M = 11000) (15) 

where a single excitation is initially present in oscillator a. 
The evolved state according to the Hamiltonian eq. il4\ may 
then be written as 

|^(t)) = a(t)|0001) + 6(t)|1000) + c(i)|0100) + d(t)|0010) 

(16) 

where, with the scaled time t = ^ and the constant a = 
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FIG. 4: Overlap between the W state with x — —y — —z = 1 /\/3 
and the evolved state of the system. The ring consists of 70 oscilla- 
tors. 

different W states one could vary the coupling constants in- 
dependently what would lead to the following Hamiltonian to 
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FIG. 5: Time evolution of the populations, |a(r)p in red, |&(t)P 
in green, and lc(r)p — jci(r)p in blue. The ring consists of 70 
oscillators. 

the system 
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In summary, we have demonstrated that it is possible to trans- 
fer quantum information with high efficiency but minimal spa- 
tial and temporal control between arbitrary sites of a trans- 
lation invariant chain of quantum systems. We have shown 
that this process works in the continuous variable regime but 
also the the single excitation regime when the system becomes 
equivalent to the dynamics exhibited by a single excitation in 
a spin chain with xy-interaction. This interaction may be gen- 
eralized to include more oscillators coupled to the chain al- 
lowing the generation of multi-particle entangled states. All 



these suggest that translation invariant chains of interacting 
quantum systems are promising candidates for the transport 
of quantum information in solid state realizations of quantum 
information processing devices. 
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